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For a system near a quantum critical point (QCP), above its lower critical dimension g!l, there is 
in general a critical line of second order phase transitions that separates the broken symmetry phase 
at finite temperatures from the disordered phase. The phase transitions along this line are governed 
by thermal critical exponents that are different from those associated with the quantum critical 
point. We point out that, if the effective dimension of the QCP, d e ff = d + z (d is the Euclidean 
dimension of the system and z the dynamic quantum critical exponent) is above its upper critical 
dimension dc, there is an intermingle of classical (thermal) and quantum critical fluctuations near 
the QCP. This is due to the breakdown of the generalized scaling relation ip = vz between the 
shift exponent if) of the critical line and the crossover exponent vz, for d + z > dc by a dangerous 
irrelevant interaction. This phenomenon has clear experimental consequences, like the suppression 
of the amplitude of classical critical fluctuations near the line of finite temperature phase transitions 
as the critical temperature is reduced approaching the QCP. 



I. INTRODUCTION 

The modern theory of critical phenomena came 
to life four decades ago, when Wilson presented 
renormalization-group tools honed to deal with divergent 
free energies pQ. Soon applied to a number of classical 
second-order phase transitions, those tools determined 
the critical exponents associated with static properties, 
such as v, associated with the divergence of the correla- 
tion length, or a, with the divergence of the specific heat, 
as well as the exponent z, associated with critical slowing 
down and necessary to describe dynamical properties. 

Later, in the mid 1970's, a pair of developments ex- 
tended the scope of the renormalization-group approach 
to quantum Hamiltonians. The first development started 
with two analyses of quantum models undergoing phase 
transitions [5J 0]. Young considered the d-dimensional 
Ising model subject to a transverse magnetic field h in 
the vicinity of the critical field h c that destroys the mag- 
netic order at zero temperature [HJ. He showed that, at 
T = 0, the model mimics the classical d + z-dimcnsional 
Hamiltonian, where z = 1 is the dynamic exponent for 
the Ising model. At nonzero temperatures, by contrast, 
the critical behavior is that of the <i-dimensional model 
with h = 0. To the same general conclusions came Hertz 
[4], who considered a broader class of models, includ- 
ing ones with z > 1. Quantum fluctuations intertwine 
the static and dynamical exponents; at the critical point, 
however, thermal fluctuations overcome them. 

The second line of development was equally important. 
By the end of 1974, its remarkable record of successes 
notwithstanding, renormalization-group theory had cov- 
ered but narrow regions of phase diagrams, and some 
feared that it would never go beyond that. Fortunately, 
such fears were unfounded. The formalism had a great 
deal more power under its hood, which became visible 
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when the Kondo Hamiltonian was diagonalized and its 
physical properties, computed [2j. The Kondo model, 
a single impurity coupled to a Fermi gas, undergoes no 
phase transition; nonetheless, when a renormalization- 
group transformation was constructed, two fixed points 
became apparent, and it was shown that the unusual 
physical properties of the model reflect the crossover of 
the Hamiltonian from the high- to the low-temperature 
fixed point [2]. 

A decade later, when the same transformation was ap- 
plied to the two-impurity Kondo problem, a QCP was 
identified at the meeting point of two competing tenden- 
cies: the coupling between each impurity and the Fermi 
gas, which tends to screen the impurity moment, and 
the antiferromagnetic interaction between the impurities, 
which tends to lock the impurity moments into a singlet 
effectively decoupled from the electron gas [6] . The fixed 
point associated with the QCP is analogous to the one 
discussed by Young [5J; in contrast with Ref. |3J however, 
Ref. E] was able to extract the phase diagram for the two- 
impurity Kondo model from the renormalization-group 
streamlines. 

The insight derived from that work offered an oppor- 
tunity to perfect the then prevalent ideas concerning the 
phase diagram of the Kondo lattice. This realized, in 
a collaboration with Japiassu and Troper the present 
author envisioned a hypothetical renormalization-group 
transformation, analogous to the one discussed in Ref. El 
argued that it should have two fixed points — a zero- 
temperature quantum fixed point and a nonzero tem- 
perature classical one — and carried out a scaling analysis 
to discuss the renormalization-group flow of the model 
Hamiltonian in their vicinity [5]- 

Dependent only on basic assumptions, the scaling anal- 
ysis proved to be as powerful as it is general. It offers an 
enlightening view of phase diagrams, one that pinpoints 
the origin of the singularities and affords easy derivation 
of scaling relations. Examples will be discussed below. 
Applied to the Kondo lattice, it showed that the quan- 
tum fixed point is unstable against thermal excitations 
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and identified special renormalization-group flow lines 
stemming out of the unstable fixed point with bound- 
aries in the phase diagram. More generally, the same 
construction provides insight and quantitative informa- 
tion on the behavior of such strongly correlated systems 
as the heavy-fermions and the high-T c compounds. Given 
the instability of the fixed point, renormalization group 
flow lines are expected to emanate from the QCP and 
carry information about it even to points that are distant 
from the quantum phase transition in the phase diagram. 

The recovery of that information is one of the excit- 
ing problems in modern condensed-matter theory. This 
paper focuses that problem in a restricted domain of the 
phase diagram, the immediate vicinity of the QCP, and 
shows that, under conditions frequently found in the lab- 
oratory, thermal fluctuations may be damped by quan- 
tum fluctuations. Both the theoretical and the experi- 
mental facets of this finding are examined. 



II. HISTORY 

In the wake of Ref. [5] numerous studies of quan- 
tum phase transitions have been reported 7-16 . Ex- 
perimental studies have focused their effects on the fi- 
nite temperature behavior of different physical quan- 
tities, sometimes even at unexpectedly high temper- 
atures. Quantum phase transitions occur in a large 
variety of systems — metals, insulators, superconductors 
[3 |H1 EH US] — whenever a parameter can be tuned by 
external means to drive a system in or out of a broken 
symmetry phase at zero temperature. Disorder can also 
be used as a control parameter [17]. And the novel area 
of cold atom systems has opened broad new perspectives, 
since the strength of the interactions, which is generally 
fixed in condensed matter materials, can now be varied 
and used to explore different types of phase diagrams 
close to a zero temperature phase transition |15j . 

Theoretically, many techniques have been used to un- 
derstand quantum phase transitions, such as, simulations 
[TO] , different types of approximations [TTJ [T51 [Jj|] , scal- 
ing theories [313] and the renormalization group, both in 
real [TU [20] and momentum space [H [TTJ [21] . These ap- 
proaches yield the phase diagrams of the relevant models, 
the behavior of physical quantities and the universality 
class of quantum critical points. In many cases, theory 
has been quite successful in explaining and even antici- 
pating experimental results [5J [7J [5] . However there are 
situations where the experimental results seem to imply 
a dimensionality for the critical quantum fluctuations of 
the system different from that expected from its crys- 
talline structure. This is the case of the ubiquitous loga- 
rithmic term in the specific heat of heavy fermions near 
an antiferromagnetic QCP. This logarithmic term arises 
in a d — 2, z — 2 universality class, but the systems 
showing this behavior in spite of anisotropic are clearly 
three dimensional |22) . 

As explained above, a distinctive feature of quantum 



phase transitions is the inextricability of static and dy- 
namics in the critical phenomena [31 H]. A clear man- 
ifestation of this entanglement between space and time 
fluctuations is expressed by the quantum hyperscaling re- 
lation [3], 2 — a = v{d + z) where a and v are standard 
critical exponents characterizing the non-analytic behav- 
ior of the free energy and the divergence of the correlation 
length, respectively [5] . The dimensionality of the system 
d appears in this relation modified by the dynamic ex- 
ponent z. It is through this relation and the scaling of 
temperature near the QCP that the dynamic exponent 
enters in the expression for the singular part of the free 
energy [23 [H] 

and ends up determining the critical behavior of static 
thermodynamic quantities. In Eq. (JTJ), h is the field con- 
jugate to the order parameter, j3 its critical exponent and 
7 that of the susceptibility. The parameter g measures 
the distance to the quantum phase transition occurring 
at g = 0. The appearance of the dynamic exponent in 
thermodynamic quantities is a unique feature of quan- 
tum phase transitions [23 25 . It turns out that if the 
effective dimension d e f f = d + z associated with a QCP 
is larger than its upper critical dimension dc, the critical 
exponents describing the quantum critical behavior are 
Gaussian, or mean field exponents [JJ . 



III. CLASSICAL AND QUANTUM CRITICAL 
FLUCTUATIONS 

In Fig. [IJ we show a schematic phase diagram of a sys- 
tem exhibiting quantum and thermal phase transitions. 
This phase diagram is typical of quantum systems above 
their lower critical dimension where the broken symme- 
try phase exists at finite temperatures. This can be rec- 
ognized, for example, as the phase diagram of the two 
dimensional Ising model in a transverse field resulting 
from a real space renormalization group approach |20j . 

Classical and quantum critical phenomena are better 
described in the language of the renormalization group, 
which associates critical points with unstable fixed points 
of scaling transformations. The phase diagram of Fig. [TJ 
contains a fully unstable zero temperature fixed point, 
the QCP that controls quantum criticality, and a semi- 
unstable fixed point at finite temperatures, the thermal 
critical point (TCP) governing the classical critical be- 
havior along the line of finite temperature phase transi- 
tions, Tc(g). The arrows in Fig. [I] show the flow of the 
renormalization-group equations. Notice that the RG 
flow along the critical line is away from the QCP and 
towards the TCP. 

In the figure, g is a control parameter such that g = 
at the QCP. In the specific case of the transverse Ising 
model, g — h — he, with he = [H/ J)c the critical ratio 
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of the transverse field and the Ising interaction. The clas- 
sical critical exponents are obtained from an expansion 
of the renormalization-group equations in the neighbor- 
hood of the thermal fixed point at TCP. Since the RG 
flow along the critical line runs towards it, the TCP de- 
termines the critical behavior along the entire line. The 
QCP controls only the zero-temperature phase transi- 
tion. However, as pointed out before, quantum fluctu- 
ations affect the finite temperature physical properties 
making the properties of the QCP experimentally acces- 
sible, a point that will be further discussed below [El US] • 
The phase diagram of Fig. [T] shows clearly a gen- 
eral feature of critical phenomena in quantum systems 
above their lower critical dimension dh- Two universal- 
ity classes coexist in the problem, related to the thermal 
and quantum phase transitions governed by the TCP and 
QCP, respectively. Expansions of the RG equations close 
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FIG. 1: (Color online) Generic phase diagram for a quantum 
system above its lower critical dimension. As a function of 
the control parameter g there is a quantum phase transition 
at the quantum critical point (QCP). There is also a critical 
line Tc(g) of finite temperature phase transitions which is 
governed by a thermal critical point (TCP). The arrows show 
the flow of the renormalization group equations. 

to these fixed points yield different sets of critical expo- 
nents. In many cases of interest, the dimensionality d is 
smaller than the upper critical dimensionality dc, while 
d + z > dc; in such cases the exponents of the QCP are 
mean-field or Gaussian, while those describing the singu- 
larities along the critical line Tc (<?) are the Wilson expo- 
nents associated with the TCP, which are much harder 
to determine. 



IV. AMPLITUDE RELATIONS CLOSE TO A 
QUANTUM CRITICAL POINT 

For the sake of clarity we choose to present the ar- 
guments in this section in the context of a real system 
showing both quantum and thermal phase transitions. 
The insulating antiferromagnetic material know as DTN 
[57] , with formula NiCl 2 -4SC(NH 2 )2 presents two quan- 



tum phase transitions as a function of an external mag- 
netic field, which have been shown to be due to a Bose- 
Einstcin condensation of magnons. Also, it has a line of 
thermal phase transitions separating the broken symme- 
try planar antiferromagnetic phase from the disordered 
paramagnetic phase. 

A schematic temperature versus magnetic field phase 
diagram for this system is shown in Fig. [2j There are 
two QCP's at Hq\ and Hqi which are in the universal- 
ity class of the density-driven three-dimensional Bose- 
Einstein condensation, with dynamic exponent z = 2 
|27) . On the other hand, the thermal phase transitions 
along the critical frontier separating the planar antiferro- 
magnetic from the paramagnetic state are in the univer- 
sality class of the classical three-dimensional XY model 
|28j . This is associated with a hypothetical finite tem- 
perature three-dimensional AYfixed point also shown in 
the figure together with the expected RG flows. 
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FIG. 2: (Color online) Schematic phase diagram of the an- 
tiferromagnetic material NiCl2-4SC(NH2)2- There are two 
QCP's at Hci and Hcz which are in the universality class of 
the zero temperature density-driven three-dimensional Bose- 
Einstein condensation with dynamic exponent z = 2. The 
finite temperature transitions along the critical line belong to 
the universality class of the classical three-dimensional XY 
model. The dashed line is one of the quantum critical trajec- 
tories (QCT). 

Consider the zero temperature Bosc-Einstein conden- 
sation phase transition at the critical external magnetic 
field Hqx. A similar analysis can be carried out near 
Hc2- For small temperatures T, sufficiently close to the 
QCP, the critical line is given by H C x(T) = H C i(l + 
uT 1 /*) where u is a constant related to the magnon- 
magnon interaction |29| . and ip the shift exponent. Al- 
ternatively, the equation for this line close to Hex can be 
written as Tc(H) oc \g\^, where g = H — Hex- 

For d+z > 4 the interaction u is a dangerous irrelevant 
quartic interaction 30 acting on the Gaussian fixed point 
describing the QCP. It is irrelevant in the RG sense since 
it scales to zero close to this fixed point. It is dangerously 
irrelevant since it appears in the denominator of certain 
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physical quantities, as shown below. It is responsible 
in this case for breaking down the generalized scaling 
relation [3T] ip — vz. 

Figure [2] shows that the flow of the RG equations at fi- 
nite temperatures runs away from the QCP and towards 
the semi-stable three-dimensional XY fixed point (TCP). 
In the language of the RG, temperature is a relevant field 
which scales away from the zero temperature QCP. The 
flow is towards the semi-unstable TCP, which determines 
the thermal critical exponents. The latter, which de- 
scribe the thermal phase transitions at Hc{T) and are 
generally different from those associated with the QCP's, 
will be denoted by a tilde to distinguish them from the 
critical exponents associated with the QCP. We will thus 
refer to the quantum exponents as {a, /?, 7, v , z, • • • } 
and to the thermal critical exponents as {a, /3, 7, v, ■ ■ ■}■ 

In the first set we emphasize the special role of the 
dynamic quantum critical exponent z, which appears in 
the expression for the free energy, Eq. ([I]). The classical 
counterpart z has no special significance, since it does 
not affect the critical behavior of static thermodynamic 
quantities. 

In each set, the critical exponents are bound by scaling 
relations, for example, a + 2(3 + j = 2 and a + 2(3 + j = 2, 
which reduce the number of independent critical expo- 
nents. Particularly important is the hyperscaling rela- 
tion involving the dimensionality d of the system. In the 
classical case the hyperscaling relation is 8J 2 — a = vd. 
The quantum hyperscaling relation in turn involves 
the dynamic quantum critical exponent and is given by 
2 — a = v{d + z). The appearance of an effective di- 
mension d e tf = d + z, has important consequences, for it 
decreases the upper Euclidean dimension dc above which 
the system is described by Gaussian or mean field expo- 
nents [SJ. In DTN, since the quantum phase transitions 
at Hex and Hci have d e // = d + z = 5 > dc — 4, both 
are described by Gaussian or mean field exponents |30j . 

The above discussion appears to suggest that classical 
and quantum phase transitions in general have separate 
descriptions not interfering with each other. That this is 
not the case will be shown below, on the basis of a gen- 
eralized scaling approach [TS] (GSA) extended for quan- 
tum phase transitions [26]. More specifically, we will see 
that, in the region of small nonzero temperatures in the 
neighborhood of the QCP, if the system is above its lower 
critical dimension dh and presents a line of second order 
thermal phase transitions, the description of the quan- 
tum and thermal critical behavior may involve, in each 
case, both quantum and thermal exponents. 

In general this occurs if the quantum phase transition 
is above its upper critical dimension, i.e., if d + z > dc 
and the generalized scaling relation ip = vz is violated 
due to a dangerous irrelevant interaction [21 j . which for 
the BEC problem considered here is the magnon-magnon 
interaction 1 30 . In this case, both quantum and thermal 
exponents are necessary to describe the physical behavior 
along the quantum critical trajectory, g — 0, T — > 0. 

For the Bose-Einstein condensation of magnons driven 



by an external magnetic field, the transverse susceptibil- 
ity x± {H, T) plays the role of the order parameter sus- 
ceptibility for this transition. Close to the QCP at Hco 
(Hco = Hci or Hco — Hci) this susceptibility has the 
scaling form [5J 

x± (F,r)«|ff(T)ro(i^p), (2) 

where g(T) = \H — Hc(T)\/Hco is the normalized dis- 
tance to the critical line g(Tc) = and Hq(T) = 
Hcofi + uT 1 /^) in the phase diagram. 

Continuity imposes the following asymptotic behaviors 
of the scaling function Q(t) (t = T/\g(T)\ uz ) [8j|3l]: 

• Q(t = 0) = constant. This guarantees that at zero 
temperature we obtain the correct expression for 
the quantum critical behavior of the transverse sus- 
ceptibility, x±(H,0) cx |g|~ 7 , such that, the order 
parameter susceptibility diverges at zero tempera- 
ture with the quantum critical exponent 7. 

• Q (t — > 00) oc t v . The exponent y — (7 — 7)/^ 
yields the correct thermal critical behavior of the 
transverse susceptibility close to the critical line 
g(T) = 0, i.e., X JH,T) cx A ± (T)\g(T)\-^ as 
T^Tc(H). 

The amplitude A±(T) cx j , (7-7)/"2 i s determined by both 
classical and quantum critical exponents. Then, close to 
the critical line, the order parameter susceptibility di- 
verges as 

X ± <*A ± (T)\g(T)r (3) 

with the correct thermal critical exponent 7. 

We now consider a special trajectory in the phase di- 
agram of Fig. [2] we let the system sit at the QCP and 
lower the temperature. Along this quantum critical tra- 
jectory (QCT), g(0) = 0, i.e., for H = H C {T = 0) = H ca 
and T — > 0, with help of the expression for the amplitude 
Aj_(T), we find from Eq. ^ that the transverse suscep- 
tibility diverges as, 

xJJJco.TjaT^uT*)-* (4) 

Then, it follows that for d + z > dc in the presence of 
a dangerous irrelevant interaction u breaking down the 
generalized scaling relation, such that ip ^ vz, both ther- 
mal and quantum critical exponents are required to give 
the correct critical behavior along the QCT. Notice how- 
ever, that for ijj = vz, the exponent 7 cancels out in 
the above expression and the divergence of the order pa- 
rameter susceptibility, x±(Hco,T) cx T ~>l vz is governed 
only by the critical exponents associated with the QCP. 
The temperature dependence of this result coincides with 
that of a purely Gaussian (u — 0) theory, in which the 
quartic interaction u is neglected (SJ. The equality be- 
tween the crossover and the shift exponents vz — ip, i.e., 
the generalized scaling relation is expected to hold for 
d + z < dc- 
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For the correlation length, a similar analysis yields 

Z<xA L (T)\g(T)r (5) 

with A L {T) = yC*-")/"*. Along the quantum critical 
trajectory, 



(6) 



and in general we find the same interference between clas- 
sical and quantum critical behavior for ip ^ vz. Again, 
for if> = vz the temperature dependence of £(i?coj T) re- 
duces to that of the purely Gaussian case, £{Hco,T) oc 
T^/ z . Notice in Eqs. Q and ^ that u appears in 
the denominator, justifying its dangerous irrelevant char- 
acter. Equations Q and ^ and their respective am- 
plitudes give rise to a new effect, namely the quantum 
suppression of classical fluctuations, if the thermal expo- 
nents are non-Gaussian. 

The specific heat requires a more careful analysis since 
it vanishes in the zero temperature axis. Its scaling be- 
havior is obtained from the scaling expression for the 
singular part of the free energy close to the critical line 
which is given by, 



f S KT^\g(T)f 



(7) 



The most singular contribution for the specific heat close 
to this line is given by, 



C/TkT^\ 9 {T)\- & (g'(T)Y 



(8) 



where g'{T) = dg/dT = -{uHcq/^T 1 ^- 1 . Along the 
QCT, if ip = vz, the temperature dependence of this 
contribution coincides with that of the purely Gaussian 
result ( u = 0) [8i, namely C/T oc T < - d ~ z ^ z . 

We have used in this section a generalized scaling ap- 
proach to show that for ip ^ vz, the critical behavior of 
physical quantities in the neighborhood of a QCP from 
which emanates a line of second order thermal phase 
transitions depends on both quantum and thermal expo- 
nents. Even along the quantum critical trajectory, g = 0, 
T — > 0, the classical and quantum critical behaviors are 
intermingled, provided only that ip ^ vz. 



V. GAUSSIAN APPROXIMATIONS FOR THE 
CLASSICAL CRITICAL BEHAVIOR 

The full problem schematized in Fig. [T] is of very diffi- 
cult solution. The real space RG in general gives global 
phase diagrams and all the relevant fixed points in the 
problem [20] . However, in this approach, it is hard to ob- 
tain controlled expansions that yield the correct critical 
exponents associated with these different fixed points. 
And for electronic systems the real space RG has only 
very limited success [8]. Momentum space RG [21] and 
self-consistent approximations [18j on the other hand give 
a better treatment of strongly correlated electronic sys- 
tems. The former |21j describes correctly the unstable 



fixed point associated with the QCP; unfortunately it is 
limited to the neighborhood of the fixed point. Although 
it describes the correct temperature runaway flow away 
from the QCP, it fails to yield the thermal Wilson fixed 
point to which the flow along the critical line is attracted. 
The self-consistent approach [18j is a Gaussian theory 
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FIG. 3: (Color online) The order parameter susceptibility 
close to the critical line Hc(T) for fixed temperatures as- 
suming that the thermal exponents are those of the three- 
dimensional XY model. The curves have increasing temper- 
ature from left to right and extend to (H — Hc{T)) /Hco = 
5 x 10~ 2 . See Fig. Blfor comparison. 
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FIG. 4: (Color online) The order parameter susceptibility 
close to the critical line Hc(T) for fixed temperatures as- 
suming Gaussian exponents for the thermal phase transitions. 
The curves have increasing temperature from left to right and 
extend to (H - H C (T))/H C0 = 5 x 1CT 2 . 

which yields the correct quantum critical exponents if the 
QCP effective dimension is larger than the upper critical 
dimension. However, it is inadequate to treat the thermal 
phase transitions and yields Gaussian thermal exponents. 
In contrast with the limitations of other approaches, the 
GSA emerges as the perfect tool to understand the sub- 
tleties of quantum phase transitions. It is particularly 
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useful to explore the intermingle of classical and quan- 
tum fluctuations near a QCP. 

Both the momentum space RG and the self-consistent 
approach to quantum phase transitions yield the shift ex- 
ponent ip — z/(d + z — 2) , showing clearly the breakdown 
of the scaling relation ip — vz for d + z > dc- The mo- 
mentum space RG shows that this breakdown is due to 
a dangerous irrelevant quartic interaction acting on the 
Gaussian QCP. In the case of DTN, this is the magnon- 
magnon interaction u. 

In this section we show that the results of the mo- 
mentum space RG and of the self-consistent approach 
for the critical behavior along the quantum critical tra- 
jectory are recovered when we assume that the thermal 
exponents take thermal Gaussian values. 

We are interested in the case d + z > dc, i- e., in 
Gaussian quantum critical exponents. Consider Eqs. Q 
and ([6| for the order parameter susceptibility and cor- 
relation length, respectively. Along the quantum critical 
trajectory, taking 7 = 7 and v = v we get for the order 
parameter susceptibility, 



X ± oc (uTi-y 
and for the correlation length, 



(9) 



(10) 



which are the results from the momentum space RG and 
self-consistent approach. These results are clearly dis- 
tinct from those given in Eqs. Q and ([6|, respectively. 

More importantly, the present analysis shows that sub- 
stitution of Gaussian values for the thermal (tilde) crit- 
ical exponents leads to temperature independent ampli- 
tudes, since the quantum exponents are Gaussian. This 
follows immediately from Eqs. ^ and ^ and the cor- 
responding amplitudes. Thus, if the thermal exponents 
turn out to be Gaussian, there will be no suppression of 
classical fluctuations in the vicinity of the QCP. This is 
illustrated in Figs. [3] and [4] 

The analysis for the specific heat is much more subtle. 
Differently from other critical exponents, the quantum 
and thermal Gaussian exponents a and a are different 
and determined by the corresponding hyperscaling rela- 
tions. For example, they yield (a — a)jvz = 1. We find 
that the most singular contribution to the specific heat 
arising from the free energy, Eq. ([7]), assuming that the 
tilde exponents are Gaussian thermal exponents is given 
by, 



C/Tcx^T" 1 !^! 1 )!- 



(11) 



Notice that in this case classical fluctuations are sup- 
pressed even with Gaussian thermal exponents (ip < 2). 
Along the QCT, using a = 1/2 in d = 3, we find that 



C/T oc u 3/2 T^ 



(12) 



in this case in the numerator. Eq. ( 12 ) is the expected 
result for the specific heat along the QCT for d = 3 as- 
suming that thermal exponents are Gaussian and taking 
into account the dangerous irrelevant interaction u. How- 
ever, the purely Gaussian result [5], which neglects the 
quartic interaction u ( u = 0), namely, C/T cx T^ d ~ z ^ z 
is in general more singular as T — > then the contribu- 
tion above, Eq. (12). For the specific heat, whether we 



consider the purely Gaussian result, which neglects the 
interaction u, or Eq. (12 1, which takes u into account, 



we find that the most singular contribution determines 
the physical behavior along the QCT. This follows from 
Eqs. (11) and (12), the right hand sides of which van- 



Differently from Eqs. ^ and ( 10 ) for the susceptibility 
and correlation length, the quartic interaction u appears 



ish for u = and can be safely neglected with respect 
to the purely Gaussian (u = 0) contribution for the spe- 
cific heat along the QCT. This is not the case for the 
susceptibility and the correlation length, which diverge 
for u — 0. Finally, we point out that the same tempera- 
ture dependence of the pure Gaussian theory (u = 0) for 
C/T is obtained from Eq. ( 12 ), if we let ip — vz on the 
right-hand side. 



VI. BOSE-EINSTEIN CONDENSATION IN 
MAGNETIC SYSTEMS 

One consequence of the difference between thermal and 
quantum critical exponents, as we have seen, is the sup- 
pression of the amplitude of critical fluctuations along 
the critical line as the critical temperature is reduced in 
the vicinity of a QCP. Besides, even along the quantum 
critical trajectory, the thermal exponents appear in the 
expressions for the critical behavior of different quanti- 
ties. These are strong results with clear experimental 
consequences. In this section we discuss the implication 
of these results for magnetic systems presenting a Bose- 
Einstein condensation of magnons as a function of the 
magnetic field. 

This is the case of DTN, which has a phase diagram 
[27] similar to that shown schematically in Fig. [2] As 
pointed out before, at Hq\ and Hc2, the quantum phase 
transitions of this three dimensional system are in the 
universality class of the density-driven Bose-Einstein con- 
densation, with dynamic exponent z = 2. Since the ef- 
fective dimensions of the QCP's at these critical fields 
are d e f / = 3 + 2 = 5 which is larger than dc = 4, the 
quantum phase transitions are described by Gaussian or 
mean-field exponents. In particular v = 1/2, such that, 
vz = 1. 

The shift exponent of the critical lines emanating from 
the QCP's is given by ip = z/(d + z — 2) = 2/3, as re- 
sults from the momentum-space RG analysis in Ref. 1211 
This is clearly distinct from the crossover exponent vz — 
z/2 = 1. The dangerously irrelevant magnon-magnon in- 
teraction u is responsible for the breakdown of the scaling 
relation ip = vz in this case }3Q] . The exponents we have 
just obtained describe well the critical behavior; in par- 
ticular, -0 = 2/3 fits very well the critical lines close to 
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the critical fields Hci and Hci 

By contrast, the thermal or tilde exponents describing 
the thermal phase transitions along the finite tempera- 
ture critical line are those of the three-dimensional XY 
model. These Wilson critical exponents are well known 
[5H[55]. They have been obtained by different methods, 
from an RG e expansion or numerical methods [34] . and 
verified experimentally [35 . These thermal exponents 
are clearly distinguishable from the Gaussian ones asso- 
ciated with the QCP's. The specific heat shows a clear 
A anomaly at the critical line associated with a critical 
exponent a 0. We therefore have an opportunity to 
test the results in the previous sections. 

Taking into account that the thermal critical expo- 
nents are different from the Gaussian exponents associ- 
ated with the QCP's, we obtain a divergent order param- 
eter susceptibility along the quantum critical trajectory, 
X ± {Hco, T) oc u-WT-^+T"/ 2 ), where we used the Gaus- 
sian value for the susceptibility exponent of the QCP, 
7 = f . Since 7 « 1.32 |34j £35] for the three-dimensional 
XY model, this yields a contribution more singular than 
that of the momentum space RG and self-consistent ap- 
proach, namely xAH C q,T) cx (wT) - ^ = (11T)- 3 / 2 . 

For the specific heat we find, C/T cx u 2 T 3 / 2 logT 
along the quantum critical trajectory, g(T = 0) = 0, 
T —> 0. We used a ~ for the specific heat exponent 
of the three-dimensional XY model [5H [55] , which im- 
plies a logarithmic dependence in this case. However, 
this term is less singular then the purely Gaussian, or- 
der zero in u contribution to the specific heat, namely, 
C/T cx T( d ~ 2 )/ 2 = T 1 / 2 . 

Let us now investigate the behavior of the thermody- 
namic quantities at the critical line for sufficiently low 
temperatures near the QCP's. The influence of the quan- 
tum exponents due to the proximity of the quantum 
phase transition is contained in the expressions for the 
amplitude. Since the quantum critical exponents are 
Gaussian, we obtain for the thermal critical behavior of 
the order parameter susceptibility, 



and for the specific heat, 



H-H C (T) 



H, 



C/T cx u 2 T 



■;+a 



CO 



H - H C {T) 
Hco 



(13) 



(14) 



where for the three-dimensional XY model, 7 w 1.32 and 
a f=a implying a lambda type anomaly for the specific 
heat [SI [35] along the critical line. For completeness, we 
quote other critical exponents for the three-dimensional 
XY model, u 0.67 and f3 « 0.34 [51 [55]. 

For the analysis of the experimental data on the critical 
behavior along the critical line we can adopt the following 
procedures. Suppose that, to approach the critical line 
in Fig. [2j we decrease the temperature for a fixed mag- 
netic field Hi > Hco- We then stop our measurements, 



say of the order parameter susceptibility, at a given arbi- 
trary temperature T\ , so distant from the thermal phase 
transition that [T ± - T c (Hi)]/H C o = S. 

Next, we repeat the measurements for another fixed 
magnetic field iJ 2 > Hi, stopping at another temper- 
ature T 2 so distant from the critical line that [T 2 — 
Tc{H2)]/Hco = 5. When this procedure is repeated 
for several magnetic fields, we find from Eq. ( 13 ) that 



the susceptibility at the temperatures Tj such that [Tj — 
T c (H t )]/H C o = Sis given by, X ±(Ti) oc lf~ 1 5 where 8 
is a fixed known constant. In this way using the values 
of x±[Ti) we can determine the exponent 7. 

Alternatively, we could increase the magnetic field at a 
fixed temperature T\. We would then stop the measure- 
ments at a field Hi so distant from the thermal phase 
transition that [H\ — Hc(Ti)}/ Hco = 8. We would re- 
peat the measurement for another fixed temperature T 2 
and stop at a distance from the critical line, such that 
[Hi - H c (T 2 )]/H C o = S. The susceptibility at the fields 
such that [Hi — Hc(Ti)]/Hco = 8 would then be given 
by x± °c T^~ 1 5 where 5 is a known constant. 

Since, in practice, the susceptibility never diverges, we 
can take the temperatures Tj as the transition temper- 
atures themselves. In this case we can also work with 
the amplitudes given directly in terms of the magnetic 
field. Since Tj = [(Hi — Hco) /uHco]^ , substitution 
on the right-hand side of Eq. (13) leads to x±(Hi) oc 

j -0(7-1) 

Fig. [3] is then the amplitude function. 

We have carried out the procedure outlined above for 
the DTN system [27] and also for Sr 3 Cr 2 8 [37] which is 
another magnetic system candidate for BE condensation 
[57] . The results for DTN are shown in Fig. [5] using the 
data near Hci from Ref. 155] The figure shows the specific 
heat at the peak temperatures, i.e., at the temperatures 
of the thermal phase transitions for different magnetic 
fields. The line through the points is C oc T 5 / 2 , as ex- 



[(Hi — Hco)/uHco) ■ The envelope of the peaks in 



pected from Eq. ( 14 ) for 5 = 



The fully Gaussian results, i.e., results obtained with 
thermal gaussian exponents, give a stronger suppression 
of the specific heat fluctuations namely, C a T 3 . This 
seems to be the case for Sr^Cr20 8 , which is shown in 
Fig. [6] here, a T 3 law fits very well the amplitude of 
the specific heat in Ref. [57J for various magnetic fields 
close to the critical field. In this case, however, the T 3 
dependence raises a doubt: it is not clear whether we are 
observing critical fluctuations or measuring the phonon 
contribution for the specific heat. 

In any case, Sr^Cr20 8 is an isotropic magnetic system 
[57], unlike DTN, which has a condensed state of the 
XY type. In the inset of Fig. [5] we show the data for the 
specific heat of DTN along the QCT at Hex- This follows 
approximately a T 3 / 2 behavior, in agreement with our 
conclusion that the purely Gaussian contribution (u = 0) 
C/T cx T( d ~ z )/ Z dominates the thermal behavior of the 
specific heat. 

Finally we discuss the critical behavior of the longi- 
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FIG. 5: (Color online) The amplitude of the specific heat 
peaks at the finite temperature phase transitions into the 
planar antiferromagnetic order or Bose condensed phase of 
DTN. The curve is a plot of C oc T 5/2 obtained using a « 
as appropriate for the three-dimensional XY model. The in- 
set shows the specific heat at the QCP at Hoi- Points are 
data from Ref. 



gitudinal susceptibility vanishes at the zero temperature 
phase transition. Mean-field would give rise to at most a 
jump in the longitudinal susceptibility at the transition. 

At T = when the external magnetic field is larger 
than the upper critical field Hc2, the system is fully 
polarized and the scaling function K(t — 0) = 0. As 
temperature increases for H > Hc2 and below the line 
T < Tq = \g\" z , this scaling function can be easily 
guessed, K(t) = exp(— Tq/T), since the system is gapped 
in this region of the phase diagram. 

For finite temperatures, close to the critical line 
Hq{T), we obtain \\\ from the second derivative of the 
free energy with respect to the magnetic field. Using 
Eq. 0, we get, 



XII 



d 2 fs 
dH 2 



(x T~ 



H - H C (T) 



Hco 



(16) 



Along the quantum critical trajectory, then, we get x\\ °^ 
T{ & - a )/ vz (uT 1 W)- & . For DTN we expect to find X \\ oc 
VTlnT along the QCT, since a « 0. 



VII. DISCUSSION 
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FIG. 6: (Color online) The amplitude of the specific 
heat peaks at the finite temperature phase transitions in 
Sr 3 Cr 2 O s . The curve is a plot of C oc T 3 with or; 1/2, as 
appropriate for the three-dimensional Gaussian model. Points 
are data from Ref. [371 



tudinal susceptibility x\\ 01 the magnetic Bose system. 
At zero temperature, Xll oc d 2 fs/dH 2 scales as the com- 
pressibility of the Bose gas, 



X\\ oc \g\ 



M" 2 



(15) 



Notice that, in three dimensions, the quantum Gaussian 
exponent a — 2— v(d+z) = —1/2 is negative and the lon- 



Although discussed in the context of insulating mag- 
netic systems presenting a Bose-Einstein condensation of 
magnons, the results in this paper should apply to any 
quantum phase transition. In particular, they also de- 
scribe quantum criticality in metallic systems. For exam- 
ple, the equation for the amplitude of the specific heat 
at a thermal phase transition near a QCP, 

holds near any superconductor quantum critical point 
with dynamic exponent z = 2 with the thermal super- 
conductor transition in the universality class of the three- 
dimensional XY model. 

The present approach relies only on universal features 
of continuous quantum phase transitions: we only as- 
sume that, near the transition, a characteristic length 
and a characteristic time contain all relevant informa- 
tion on the system. Any difference in the nature of the 
materials, whether metal or insulator, or in the phase 
transition itself, whether magnetic or superfluid, will be 
accounted for by the various critical exponents, in special 
by the dynamic exponent z, which distinguish different 
universality classes. 



VIII. CONCLUSION 

We have used a generalized scaling approach to study 
the critical behavior of physical quantities near a quan- 
tum phase transition. We considered systems above 
their lower critical dimension, such that a line of finite- 
temperature continuous phase transitions emanates from 
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the QCP. Our approach has shown that for quantum 
phase transitions described by Gaussian exponents, i.e., 
with d+z > dc, where dc is the upper critical dimension, 
classical and quantum fluctuations are intermingled near 
the QCP. The breakdown of the generalized scaling rela- 
tion ip — vz, between the shift and the crossover expo- 
nents, signals the intertwining. In general this breakdown 
is due to a dangerous irrelevant interaction on the Gaus- 
sian fixed point describing the zero-temperature phase 
transition. 

Our approach brings into the stage of quantum crit- 
icality, the thermal phase transitions themselves. We 
find that thermal critical fluctuations close to a QCP are 
suppressed by quantum effects and their study allows a 
deeper understanding of quantum criticality. As future 
experiments progress towards lower temperatures, more 
detailed comparisons with theory will become possible 
allowing finer distinction between alternative scenarios. 

We have shown that the momentum space RG and 
the self-consistent approach yield quantum-criticality re- 
sults that follow from the GSA under the assumption 
of Gaussian exponents for the thermal phase transitions. 
Whether this is a sufficiently good description of the ther- 
mal phase transitions and, as implied by our analysis, of 
quantum criticality itself can be tested experimentally. 
The alternative assumption, that the thermal exponents 
have the expected Wilson character, has stronger exper- 
imental consequences. The amplitude of critical thermal 



fluctuations will be reduced and the physical behavior 
along the quantum critical line will show deviations from 
the naively expected quantum critical exponents. 

We found the behavior of the specific heat near the 
QCP to be different from that of other physical quantities 
such as, for example, the order parameter susceptibility 
or the correlation length. The former has its amplitude 
suppressed near the QCP even assuming Gaussian ther- 
mal exponents. Besides, the corrections for this quantity 
due to the irrelevant interaction can safely be neglected 
when compared to the, in general more singular, purely 
Gaussian (u = 0) result. This is definitely not the case 
for other physical quantities where the corrections in u 
to the purely Gaussian result have u in the denominator 
due to its dangerously irrelevant character. 
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